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ABSTRACT 

We have obtained accurate CCD narrow band Hj3 and Ha photometry of Giant HII 
Regions in M 33, NGC 6822 and M 101. Comparison with previous determinations 
of emission line fluxes show large discrepancies; their probable origins are discussed. 
Combining our new photometric data with global velocity dispersion (a) derived from 
emission line widths we have reviewed the L{Hj3)) — a relation. A reanalysis of the 
properties of the GEHRs included in our sample shows that age spread and the su- 
perposition of components in multiple regions introduce a considerable spread in the 
regression. Combining the information available in the literature regarding ages of the 
associated clusters, evolutionary footprints on the interstellar medium, and kinemati- 
cal properties of the knots that build up the multiple GEHRs, we have found that a 
subsample - which we refer to as young and single GEHRs - do follow a tight relation 
in the L — a plane. 
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1 INTRODUCTION 

Accurate photometry of Giant HII Regions (GEHRs) is a 
key source of information about the energetics of these large 
star-forming regions. The luminosity emitted through the 
recombination lines of ionised hydrogen (H a , H^) allows to 
estimate the number of ionising photons which in turn es- 
tablish the requirements for the stellar population respon- 
sible for such an ionising flux. Several reviews have alread y 
discusse d the se properties, s uch a s Kennicutt (1979; K79), 
Shields (|l990j) and Melnick fll99S|), 



tional explanation for this observed kinematical behaviour. 
Other authors confirmed the existence of a correlation, al- 
though they found different slopes for the linear regression. 
New models were produced to account fo r the underlying dy- 
namics. Turbulence (L oc a 3 , Hippel ein ( 1986)) and multiple 
shells (L oc a 6 , Arsenault & Roy ( [19881)7 were introduced 
in the models. Fuentes Masip et al. (2000) have obtained 



These emission lines had been fo und to have supersonic 
widths by Smith & Weedman ( |l970| ) , ( |l97l[ ) . Such behaviour 
is puzzling, as supersonic motion cannot be susta ined with- 
out a source of energ y to fuel it. Melnick ( 1977 ) and later 
Terlevich & Melnick (1981) found that the gas velocity dis- 
persion (a) measured from a recombination line in giant HII 
regions correlates well with the total luminosity (L) emitted 
in the respective line, and also with the overall radius of the 
nebula (R). In a logarithmic plane, total flux a nd nebula 



radius follow a lin ear relation with the line width ( Terlevich 



and Melnickl981) and the following approximations hold: 
L oc a and R oc a 2 . These dependences resemble much 
those of virialised systems, such as globular clusters and 
elliptical galaxies, and lead the authors to suggest a gravita- 
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imaging spectroscopy of the irregular galaxy NGC 4449 with 
which they detected and analysed its giant HII regions. They 
studied the correlation of size and luminosity with the ve- 
locity dispersion and found a regression for the subsample 
that showed supersonic motions, high surface brightness and 
single Gaussian emission lines. This linear regression in the 
log plane had a slope of 3.8 ±1.2 which is consistent with a 
gravitational origin of the supersonic motions. 

The intriguing aspect of the L-a relation is that al- 
though we know that massive stars produce strong winds 
and in many cases we can see features that look like shells, 
the scaling laws seem to suggest that gravity may play a 
central role. Attempts to determine the main broadening 
mechanism in HII regions with supersonic linewidths were 
made for the closest and most luminous GEHRs within the 
Local Group, 30 Doradus in the Large M agella nic Cloud and 
NGC 604 in M 33. Chu & Kennicutt fll994j) used longslit 



echelle spectroscopy data on several positions of 30 Doradus 
and found that its global supersonic linewidth breaks up 
into a large number of smaller supersonic structures, which 
combine with a smooth underlying turbulent velocity field. 
They suggested that gravitation could not contribute signifi- 
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cantly to the observed nebular kinematics. This finding was 
partly confirmed by Melnick, Tenorio-Tagle and Terlevich 
(1999) who obtained a single longslit echelle spectrum with 
higher resolution and signal to noise ratio. The slit used 
to obtain the spectrum, passed through several important 
features which are analysed in their paper. However, they 
find no discrete component with supersonic velocity widths, 
beyond a very broad component (a ~ 45 km s _1 ) found to- 
wards all lines of sight. This does not seem to b e the case 
for NGC 604. Muiioz-Tunon et al. ( |l996; MTCl] ), analysed 
high spatial and spectral resolution Ha emission line data 
on NGC 604 and concluded that although there is a contri- 
bution from shells and filaments to the observed gas kine- 
matics, the overall supersonic behaviour is dominated by 
unresolved components, observed throughout the whole in- 
tensity range. By modelling the expected behaviour of the 
emission line supersonic width and intensity of a GEHR they 
suggested evolutionary effects to explain their di fferen t ob- 
served behaviour. In a similar study, Yang et al. (1996) con- 
cluded that a substantial fraction of the broadening must be 
due to gravity. 

The regression among global characteristics of GEHRs 
also holds for HII Galaxies. These are luminous star forming 
regions that dominate the galaxian energy output. The rela- 
tion between a distance independent parameter, such as the 
linewidth, and the luminosity of the object allows us to use 
them as distance indicators. The strength of the flux emit- 
ted in the recombination lines of these objects makes them 
available as far as z ~ 3, and even further wi th th e new 
generation large telescopes. Lowenthal et al. (1997) have 
already obtained Keck spectroscopy of a sam ple of z ~ 3 
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galaxies from the Hubble Deep Field project (Williams et 
al. 19960. They found that these objects seem to lie at the 
bright end of the luminosity - size relation present for HII 
galaxies although the starburst origin of t heir o bserved spec- 
tra remains to be prove n. Melnick et al. (2000), makin g use 
of Guzman et al. ( 1997) together with Koo et al. (1996) and 
Pettini et al. (199£) data for high z, have re-analysed the 
luminosity - linewidth regression for HII galaxies and found 
a good correlation up to z ~ 3, in spite of several systematic 
effects that might be increasing the observed scatter in the 
regression. 

Previous studies dealing with the L — a regression for 
the Giant HII regions have used datasets from a variety 
of sources and detectors, which may lead to large internal 
errors of the samples. The largest and more homogeneous 
available photometric data in the literature for the Ha emis- 
sion line in GEHRs is K79, w hich is used by Hippelein (198f ) 
and Arsenault & Roy (198£ ) in their st udies. H,g photometry 
was pe rformed by Me lnick (1979; M7£) and later by Melnick 
et al. (1987; MMTG). These data are discrepant, reasons 
for that could be the use of non-linear detectors, the lack of 
good standards for narrowband emission line photometry, or 
uncertainties in the extinction corrections. 

The disagreements between previous luminosity estima- 
tions of GEHRs suggests that errors in the determination of 
fluxes may be a main contributor to the observed spread in 
the L — a regression. The goal of the present photometric 
study is to achieve a standard set of Ha emission fluxes from 
nearby extragalactic giant HII regions to use as the 'local' 
reference system. To this end we have performed CCD pho- 
tometry on a set of GEHRs from galaxies with relatively well 



Figure 1. Finding chart for NGC 588. The axes show the angular 
scale, estimated from converting the pixel to arcsec scale for the 
TeK CCD at the JKT focal plane. The coordinates zero point 
coincides with the origin of the curve of growth method used to 
find the total flux of the region. 



Figure2 



Figure 2. Same as figure m for NGC 592. 



known distances, in order to calculate their absolute fluxes. 
The observations are described in Section bl the data reduc- 
tion is described in Section ^| and Section^ introduces the 
flux calibration process designed for this work. The proce- 
dures used to obtain the total fluxes in each band are listed 
in Section |E] and the results, plus comparison with previous 
studies are shown in Section |^. 

Extinction maps are presented in Section |^. The relation 
between the luminosity and width of the emission lines is 
discussed in Section H and conclusions are drawn in Section 



2 OBSERVATIONS 

CCD frames of the selected Giant HII Regions were obtained 
during seven nights at the 1.0m Jacobus Kapteyn Tele- 
scope (JKT) in the Observatorio del Roque de los Mucha- 
chos (ORM) La Palma, Spain, in August 1997. Five out 
of the seven nights were of excellent photometric quality; a 
sixth one, although not photometric, was useful too. The last 
night of the run had variable sky conditions so the data was 
not included in the analysis. The average seeing conditions 
were monitored from the image quality of short exposures 
of the flux standard stars, and an upper limit was also es- 
timated from stellar sources present in our longer exposed 
nebular frames. These estimations, and the atmospheric ex- 
tinction measurements available from the Carlsberg Merid- 
ian Circle, also located at La Palma, are listed in Table [ij 

The target objects, together with the flux standard stars 
were observed using four narrow-band filters. These are cen- 
tred in Ha (6561 A), H/3 (486lA), and their adjacent con- 
tinua (6656A and 5300A respectively). The main character- 
istics of the set of filters used are given in Table |i[ Three of 
them are from the ING bank of filters. The H/3 continuum 
filter was kindly provided by Don Pollaco, as the presence of 
nebular emission lines, such as [OIII]AA 4959,5007 A made 
the available continuum filters useless. 

The target objects from three different Local Group 
galaxies, namely NGC 6822, M 33, and M 101, were widely 
spread in right ascension. This allowed to sequentially ob- 
serve them quite close to their passage through the local 
meridian, reducing the effects of large airmass corrections. 
On the other hand, the objects were observed on different 
days at different bandwidths, which complicates the cali- 
bration. Table ^ gives the observing log and figures |l] to [To| 
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Figure 3. Same as figure [j] for NGC 595. 
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Figure4 

Figure 4. Same as figure U for NGC 604. 
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Figure 5. Same as figure m for NGC 5447. 
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Figure 6. Same as figure m for NGC 5461. 
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Figure 7. Same as figure fJ for NGC 5462. 
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Figure 8. Same as figure fJ for NGC 5471. 
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Figure 9. Same as figure [j] for NGC 6822 I. 



FigurelO 



Figure 10. Same as figure | for NGC 6822 III. 



Date 


Seeing 


Extinction 




(average) 


(mag.) 


09/08/97 


0'.'90 


.11 


10/08/97 


0'.'86 


.11 


11/08/97 


0'.'91 


n/a 


12/08/97 


l'/02 


.11 


13/08/97 


l'.'OO 


.11 


14/08/97 


Of.' 96 


.19 


15/08/97 


> 1'.'12 


.18 



Table 1. Sky conditions for the La Palma sky during the observ- 
ing run. Seeing conditions were measured from stellar sources on 
the CCD images. The extinction data is provided by the Carls- 
berg Meridian Circle (La Palma) 



Filter 


ING# 


A c 


AA 


Peak Tr. 


Ha 


60 


6561 


38 


33.9 




67 


6656 


11 


47.5 


H/3 


81 


4861 


54 


59.7 


H/3 corl t 




5300 


50 


58.5 



Table 2. Description of the filters used. Column 2 indicates 
the corresponding number in the Isaac Newton Group (ING) 
database. Column 3 lists the central wavelength of each filter 
(in A), column 4 the bandpass width (in A), and column 5 their 
peak transmission values (in percentage) 



Object 


Ha 




Ha co „t 


H/3 


H/3 CO ni 


Mpp coo 
IN KjKj Ooo 








11 


1 -1 
11 


Mpri rr.no 


1 1 
11 




1 1 


1 tr 
10 


1 tr 
10 


in Oyo 


JLU 




i n 


i n 

1U 


11 


IN DU4 


Q 




Q 

y 


lo 


lo 


IN V. i v_/ 044 / 


la 




1 Q 

lo 






NGC 5455 


15 




14 






NGC 5461 


11 




11 






NGC 5462 


11 




11 






NGC 5471 


10 




10 






NGC 68221 


9 




9 


10 13 


11 14 


NGC 6822III 


9 




9 


10 13 


11 14 


BD+28°4211 


9* 


10 11* 


9* 10 11* 


10 13 


11* 13 




12* 


13 15* 


12 13 14* 


14* 15* 


14* 15* 


BD+33°2642 


10* 


11* 12 


10* 11* 12 


10*13* 


11* 13* 




13* 


15 


13 14 


14 15 


14 15 


FeigellO 


9 




9 







Table 3. Observing log for the target and calibration objects. The 
number represents the date in which the object was observed. For 
the standard stars, a '*' indicates more than one observation in 
that day. 



show the emission line flux distribution in the sky. These 
figures include a scale bar, as a guideline for linear sizes of 
the regions analysed in this paper. 

Bias frames and evening and morning flats were taken 
every day. The importance of the double set of flats is dis- 
cussed in the following section. 



3 DATA REDUCTION 

Standard procedures using IRAF software were followed for 
the initial steps of the data reduction. Bias frames were com- 
bined and used for calculating the additive zero point cor- 
rection as usual. The correction for flat field, however, was 
more delicate as the nature of the flux calibration process 
discussed in the next section relies strongly on the flatness 
of the field illumination. 

As mentioned above, we obtained blank sky exposures 
to correct from non- uniform illumination of the CCD. A 
simple check consisting on correcting one set of flat field 
images, say the evening ones, with the ones obtained the 
morning after, showed a non-uniform image with patterns 
still present in it. Figure [n] shows a section of a flatfield im- 
age taken with the Ha filter. Two sets of annular patterns 
can be seen. Some are small with clearly defined edges while 
others are larger with difuse edges. The smaller rings are 
neatly removed but the diffuse annuli remain. Ring-like pat- 
terns are usually associated with grains of dust, the reason 
for their appearance being that they are observed severely 
out of focus; they are usually referred to as 'dust flakes'. 
The observed displacement might be explained if large rings 
are present in the filters. As the filter wheel does not park 
exactly in the same position each time the same filter is 
requested, the ring patterns are shifted. The smaller rings 
might therefore be present on other surfaces (i.e. dewar win- 
dow) The dust flakes problem has its maximum effect in the 
Ha filter, and is negligible for the other ones. 
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Figurell 



Fitting of extinction 



Figure 11. Section of the evening H a Hatfield frame. Note the 
large ring structure throughout the whole region, and the smaller, 
sharper rings. 





BD 28°4211 


BD 33° 2642 


Feige 110 


F(H a ) 


1041.30 


922.61 


304.23 


mag(H a ) 


-7.5439 


-7.4125 


-6.2080 


F(6650A) 


1070.84 


881.30 


320.41 


mag(6650A) 


-7.5743 


-7.3628 


-6.2643 


F(5300A) 


2575.10 


1901.85 


764.53 


map(5300A) 


-8.5270 


-8.1979 


-7.2085 




3484.94 


2475.61 


995.02 


mag(Hp) 


-8.8555 


-8.4842 


-7.4946 



Table 4. Calculated fluxes and magnitudes for the standard stars 
used for the calibration of the photometry. 



4 FLUX CALIBRATION 
4.1 Standard star magnitudes 

After the preliminary steps mentioned in the previous sec- 
tion, the next stage is the flux calibration of our photome- 
try. Narrowband emission line photometry is unusual. There 
is no set of standard stars calibrated as to perform a rou- 
tine reduction process. There is no analog to the UBV sys- 
tem to use with Ha, H/3 and their respective continuum 
filters. So to obtain photometry of emission line regions, like 
the GEHRs, some ad-hoc calibration has to be adopted. To 
achieve this, we have chosen three HST spectroph otometric 
stars (BD +28°4211, BD+33°2642, and FeigellO ( [Turnshek 
et al.1990)) for which flux is calibrated and tabulated every 
1 or 2 A. This high resolution allows to calculate their ab- 
solute flux (or magnitude) through every single filter of our 
system. This can be achieved integrating numerically, as the 
transmission curves of the used filters, scanned at the ORM, 
are available. We can define the absolute flux in a filter as: 



K(A ) 



/ F«(\)T(\)d\ 
J T(X)d\ 



(1) 



where Ao is the central wavelength of the filter with trans- 
mission curve T(A) and -F*(A) is the absolute flux of the 
star. 

The calibration magnitude of the observed standard 
stars for each filter was calculated as m = —2.5 logF, as 
it would be measured outside the earth's atmosphere. The 
results are listed in Table ^. 

4.2 Atmospheric extinction correction 

The disadvantage of working with these flux standard stars 
is the limited number of observations per night, as there 
is only one star per field. Nevertheless, as we had perfect 
photometric conditions during the first four nights of our 
run, we could gather the data available from these nights 
to perform the extinction correction. Also, within our own 
photometric system, the transformation from observed to 



Figure 12. Visualisation of the result of the least squares fit to 
the extinction correction. Each filter is represented by a different 
symbol: open triangles indicate the 5300 A filter, open circles the 
6650 A filter, asterisks the H/3 filter, and crosses show data for 
the Ha filter. Different stars are included in each filter group. The 
point where the extrapolation of the lines drawn meet the X = 
axis indicate the zero-point correction. Within the observational 
errors the assumption of a single slope for all filters is valid. 



calibrated magnitudes does not require the introduction of 
colour corrections. The transformation equations are written 
as 



m ca iib = m hs + Z\ + kX 



(2) 



where m ca ;,;j, is the star calibrated magnitude and m (, s is the 
observed magnitude, Z\ is the zero-point correction for each 
filter (to account for slight differences between the tabulated 
and real transmission of the filter) , X is the normalised air- 
mass (X — 1 at zenith) and k is the extinction in magnitudes 
due to the airmass X. 

The magnitudes of the standard stars were obtained 
performing aperture photometry on the observed stars, us- 
ing IRAF's version of DAOPHOT. We considered a different 
value of Z for each filter, but assumed the same value of k 
for all filters. As we are working within a small range of 
wavelengths (4860 - 6650 A) this is a reasonable approach. 
The set is solved numerically by least squares. The results 
are shown in Figure [l^ and listed in Table |E| 

The set of defined extinction corrected magnitudes 
is applied to the following nights, when the sky conditions 
are expected to be slightly different. For these nights we can 
apply the results of the calibration, where the k parameter 
is known, and find the new value of Z\, which might dif- 
fer from the previous one due to atmospheric changes. The 
values of AZ\ were checked, in order to see if conditions 
also changed during the night, which was the case for the 
last night of observations. The number of standard stars ob- 
servations were not enough to quantify this variability and 
account for it, so we decided not to use the data from the 
last night. 
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Parameter 


Fit Value 


Error 




-0.0766 


0.036 




0.2421 


0.037 


z Hf3 


0.0279 


0.046 


•^5300 


0.8166 


0.042 


k 


-0.0900 


0.026 



Table 5. Values from the least squares fit to solve for the extinc- 
tion correction. Z\ indicates the zero-point correction and k is 
the extinction per airmass unit. All quantities are in magnitudes. 



Region 


[Nil] A6548 


[Nil] A6584 


k 


Ref. 


NGC 588 


0.012 


0.034 


0.975 


1 


NGC 592 


0.045 


0.125 


0.912 


* 


NGC 595 


0.047 


0.134 


0.903 


2 


NGC 604 


0.043 


0.117 


0.907 


1 


NGC 5447 


0.023 


0.070 


0.961 


3 


NGC 5461 


0.037 


0.110 


0.941 


4 


NGC 5462 


0.040 


0.119 


0.937 


3 


NGC 5471 


0.008 


0.026 


0.987 


4 


NGC 68221 




0.018 


0.992 


5 


NGC 6822III 




0.019 


0.992 


5 



4.3 Calibrating the frames 

Having flatfielded the frames and defined a relation between 
absolute flux (or magnitude) and the number of counts de- 
tected, we can use the flux-to-counts relation established for 
the standard stars throughout the whole field of the CCD 
and apply it to our program objects. Before this, images were 
aligned, cosmic ray cleaned and had their continuum con- 
tribution subtracted following stan dard procedures. More 
details can be found in Bosch ( 1999 ). 

The flux standard stars have almost featureless con- 
tinua. Hence the flux remains almost constant throughout 
the whole filter range. This is not the case for the line flux 
emitted by the nebulae, as the emission is concentrated in 
less than one A. The factor by which we are underestimat- 
ing the line flux when assuming it has a flat distribution is 
calculated. The true flux can be obtained making use of spec- 
troscopy already published in the literature from which we 
know the central wavelength and dispersion of the emission 
line for our objects. Therefore we can construct a Gaussian 
emission line 



27TCT 



(A-Aof 
ex P( ^r~) 



2a 2 



(3) 



where Ao and a are selected according to the known radial 
velocities and observed line widths. Galaxies in the Local 
Group do not have large radial velocities, so the wavelength 
of the emission lines are close to the rest value. The typical 
width (a) of the emission lines comes from their supersonic 
velocity dispersion, ~ 30 km s _1 , which corresponds to a 
width of approximately 0.67 A for the Ha line. 

The Gaussian emission line defined in Eq.^ has unit flux 
by definition, so we can use a similar procedure to the one 
used for the continuum. In this case we have to calculate the 
throughput of this artificial line through the filter, 



Fneb(Xo) = 



J F neb (X)T{\)dX 
J T(X)dX ' 



(4) 



which yields the conversion factor, (F ne i,(Xo))~ , necessary 
to correct from the observed flux into true flux for the rele- 
vant emission line. 

The narrow Ha emission line, is not the only contri- 
bution to the measured flux in the Ha filter, due to the 
presence of the forbidden [Nil] AA6548, 6584 A lines within 
the range of the filter. In order to estimate - and correct 
for - the intensity of these lines, we have gathered informa- 
tion available in the literature about their relative intensities 
with respect to Ha and introduced in the simulation de- 
scribed above, Gaussian profiles with their respective inten- 



Table 6. Intensities for the [Nil] AA6548, 6584 A lines, relative 
to Ha, and the correcti on fac tor derived. Reference value s are 
from: ( I) Vi' lchcz et al. ( 1988 ) (21 Kwitter & Aller (|l98l|) , (3) 
Smith (|l975|), (4) Rayo et al. (|l982[), (5) Pagel et al. (p80|, (*) 



Averaged from the intensities of NGC 595 & NGC 604 (see text). 



sities. The net expected true intensity from the Ha emission 
has then been calculated. 

The line intensities from the literature are listed in Ta- 
ble |. In the case of NGC 5462 and NGC 54 47 on ly the total 
value for both lines was available in Smith (1975). To follow 
the procedure described above, we need to introduce inten- 
sities for each [Nil] line in the simulated spectrum in order 
to calculate their relative effect due to the observed redshift. 
The intensity for each of the [Nil] lines was obtained from 
the total one, considering average intensity ratios between 
both [Nil] lines (A6548/A6584 ~ 3). The value for NGC 592 
(which doesn't exist in the literature) was estimated consid- 
ering the fact that there is a rad ial gradient in the abundance 
of nitrogen detected in M 33 (Vilchez et al.1985) and that 
the region is 3.4 kpc away from the centre of the galaxy, 
quite close to the radial distance of NGC 604 and NGC 595. 
An average of the intensities of the [Nil] lines for these two 
regions was used to estimate the intensity of the lines in 
NGC 592. 



5 Ha AND H/3 PHOTOMETRY 

The photometric analysis was done with the PHOTOM, 
GAIA and ESP packages from STARLINK. These allow to 
measure flux within increasing radius from the centre and 
build a curve of growth in order to obtain the total flux of 
the GEHR in each emission line. 

The centre point for calculating the curve of growth was 
selected to coincide with the position of the ionising stellar 
cluster for the rounded, symmetrical regions. The position 
of each cluster was determined from the observations in the 
previously aligned continuum bands. In the case of multiple 
regions, the curve of growth starting point was found by cen- 
tring the global emission of the whole region. Although this 
lead to centres which did not correspond to any observed 
feature in particular, the shape of the curve of growth was 
better defined and allowed a better estimation of its asymp- 
totic value. To do this we compared the flux increase on 
three succesive rings to the estimated flux noise. When this 
increment was found to be smaller than the uncertainties, 
the aperture integration was stopped and the final flux value 
recorded. 
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Figure 13. Radial and cumulative flux density for the GEHRs in 
M 33. X-axis indicate radial distance to the region centre. Radial 
flux density profiles (solid lines) are labelled in the left y-axis, 
The incremental total flux profile (dashed lines) are labelled in 
the right y-axis. Each region is identified on the cumulative flux 
plot. The top to bottom order holds for the radial density plots. 
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Figure 14. Radial and cumulative flux density for the GEHRs 
in M 101. Axes and labels are as in figure [131 



By visual inspection of the images, the residuals due 
to over or under subtraction of stars were detected, marked 
as 'bad pixel' regions and neglected when calculating the 
observed flux. 



Figure 15. Radial and cumulative flux density for the GEHRs 
in NGC 6822. Axes and labels as in figure 



5.1 Estimation of the photometric errors 

The curve of growth procedure described previously yields 
the total flux detected for each region. The error bars shown 
in the corresponding plots ( |l3[ represent the uncer- 

tainty due to the photon counting process, determined by 
the AUTOPHOTOM package. As we are dealing with ex- 
tended objects, the errors in the determination of the back- 
ground emission are very important. To account for this, we 
estimate the noise in the photometry from the dispersion 
measured in the regions used to perform the background 
subtraction instead of using the nominal readout noise of 
the CCD. This ensures a more realistic estimate of the un- 
certainties involved in the process. 

As for the flux calibration errors, there is an unavoid- 
able uncertainty in the zero-point correction (Z\) and the 
extinction per airmass unit (k). The overall uncertainty can 
be calculated then differentiating the expression that links 
observed flux (F) and counts (N) in the detector: 

OF OF dF 

AF -dN AN+ 0Z AZ+ 0k Ak 
which leads to 



(•») 



AF 



Z + kX 
10T5- 



z+kx ln(10) 



2.5 



AZ 



Z + kX 



ln(10) 
2.5 



X Ak. 



(6) 



Dividing this uncertainty by the total flux measured, 
we can get the relative error in the flux determination: 

AF _ AiV ln(10) 



— = — + ^ rT ±(AZ + AkX) 



2.5 



(7) 



From the curves of growth (figures to ^) it can be 
seen that the relative error due to the photon counting and 
background noise is about 4%, or less. The goodness of the 
extinction correction was also shown in the parameters of 
the fit to the standard stars included in table ^ From this 
table an average error of 0.04 can be estimated for the zero 
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point correction and of 0.03 for the uncertainty in the deter- 
mination of the extinction per airmass unit. As the airmass 
of the observed objects was always smaller than 2, a rough 
average of 1.5 for X leads to an estimate for the second term 
of equation ^ of 7%. Finally a 2% error should be accounted 
for the 'dust flakes' problem in the flat fields. Hence, we 
consider that a value of 13% correctly represents - conser- 
vatively - the relative error in the fluxes determined in this 
work. 



5.2 Interstellar Extinction Correction 

Both H a and Hp fluxes were obtained for six of the observed 
regions allowing an extinction correction. We will summarise 
the standard procedure now, specifying the nomenclature 
used to avoid confusion: 

The observed flux F at a wavelength A of an intrinsic 
flux Fo after passing through a medium of optical depth r\ 
is 



Fx = F X0 e 



(8) 



The ratio of observed to true intensities at different 
wavelengths is then 

F Xl = Fx l0 
F\„ F\ 2 o 



(9) 



It is assumed that the interstellar extinction has the 
same wavelength dependence, so 

F\i_ _ F\iO e - c [/(Ai)-/(A 2 )] 

F\ 2 F\ 2 o 



(10) 



F H 



Using logarithms in base 10, and for the case of H a and 
l -c[f{H a )-f(H ff )] 



Hp we have 
F Ha F Ha o 







F H/3 o 



where C = 0.434 c, and taking logarithms 



log 



Fh, 
Fh 



(11) 



f(H )), (12) 



allows us to calculate the extinction factor C. The measured 
ratio of observed fluxes can be compared with the theoret- 
ical ratio of the same pair of lines for an HII region. Case 
B recombination theory for an ionisation bound nebula at 
10 4 K, predicts th e logarithm of the ratio of Hp and H a 
fluxes to be 0.454 (|Brocklchurstl97l|) . The / (A) values are 
obtained from Osterbrock (1989). This value of C is found 
in the literature as C(Hp), the logarithmic Hp absorption 
coefficient. 

The logarithmic H a absorption coefficient, which we 
need, can be calculated from C(Hp) as 

C( J F/ A ) = log(F A 0)-log(F A ) (13) 

so, we can rearrange Eq,|l2] 

log (F H J - log (F Ha0 ) = -C(Hp) + C(Hp) /, (14) 
where / is / (H a ) — f (Hp). Replacing Eq.[l3] for H a : 
C(H a ) = C(Hp) (1 - /) (15) 



GEHR id. 




F H p 




NGC 588 


-11.48 


-11.87 


0.00 


NGC 592 


-11.57 






NGC 595 


-10.97 


-11.44 


0.01 


NGC 604 


-10.47 


-10.93 


0.00 


NGC 5447 


-11.33 






NGC 5461 


-11.19 






NGC 5462 


-11.37 






NGC 5471 


-11.21 






NGC 6822 I 


-11.36 


-11.86 


0.03 


NGC 6822 III 


-11.32 


-11.87 


0.06 



Table 7. Results from the photometric analysis of the observed 
giant HII regions. For each region, we have listed the Ha flux in 
column 2, the H/3 flux in column 3, and the logarithmic absorp- 
tion coefficient in column 4. 



6 RESULTS 

The total flux for each region was derived from the asymp- 
totic value determined by the curve of growth method. In 
order to analyse the effects of residuals in the sky and contin- 
uum subtraction, we have performed the photometry using 
different methods for sky estimation and masking of the fea- 
tures not perfectly removed by the continuum subtraction. 

The measured Ha and Hp fluxes in the available band- 
widths are listed in Table together with our determina- 
tions for the logarithmic H a extinction coefficient. 'Nil' val- 
ues correspond to those cases in which the measured ratio 
was equal to, or slightly lower than, the theoretical one. This 
might be suggesting there is some systematic problem be- 
tween fluxes determined for both sets of filters, and we will 
discuss this in the following sections. 



6.1 Comparison with previous photometric 
studies 

The largest collection of total fluxes of GEHRs in the Local 
Group in the Hf3 bandwidth have been obtained by M79 
and MMTG. These are summarised in Table ^. We have in- 
cluded in column 2, the expected flux of Yip from K79 H a 
photometry, calculated using the theoretical Ha/Hj3 ratio 
for Case-2 recombination (Osterbrock 1989) . The discrepan- 
cies among the measured fluxes are easily noticed. There are 
several possible sources for these discrepancies. Firstly, the 
use of non-linear detectors can fail to register properly the 
large dynamical range of fluxes in GEHRs from the weak 
diffuse emission present over a large area to the strong flux 
emitted by conspicuous 'knots'. Secondly, the flux calibra- 
tion process for narrowband emission line photometry lacks 
of standard objects and a proper reference system, such as 
Johnson & Morgan ( 1953 ) for the UBV photometry. Thirdly, 
the extinction correction also varies among different stud- 
ies, and can account for some of the discrepancies in the net 
fluxes calculated. 

The comparison of our results with previously pub- 
lished photometry on these HII regions is discussed sepa- 
rately in the following four sub-sections, in order to anal- 
yse the possible sources of disagreement. The differences be- 
tween published fluxes and the ones obtained in the present 
study are presented in table H. These are calculated as 
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GEHR id. Kennicutt(1979) Melnick(1979) MMTG 

FHa F^jp F H p 



NGC 588 


-11.63 


-12.08 


-12.09 


-11.73 


NGC 592 


-11.82 


-12.27 




-11.61 


NGC 595 


-11.06 


-11.51 


-11.74 


-11.58 


NGC 604 


-10.60 


-10.95 


-11.15 


-11.17 


NGC 5447 


-11.80 


-12.25 




-12.57 


NGC 5461 


-11.50 


-11.95 


-11.89 


-12.29 


NGC 5462 


-11.82 


-12.27 


-12.37 


-12.21 


NGC 5471 


-11.45 


-11.90 


-11.97 


-12.10 


NGC 6822 I 


-11.30 


-11.75 


-12.09 


-12.48 


NGC 6822 III 


-11.23 


-11.68 


-12.05 


-12.49 



Table 8. Published photometry for Giant HII Regions in the Lo- 
cal Group. This subset includes the regions for which we have 
obtained new photometric data. The H/3 fluxes for K79 are es- 
timated from the Ho flux assuming a theoretical ratio between 
these two lines of 2.82, only for comparison purposes 



GEHR id. 


Ajf79 


Am 79 


Ammtg 




Fhc 


Fhp 


Fh/3 


NGC 588 


0.17 


0.26 


-0.10 


NGC 592 


0.22 






NGC 595 


0.08 


0.37 


0.21 


NGC 604 


0.09 


0.18 


0.20 


NGC 5447 


0.43 






NGC 5461 


0.30 






NGC 5462 


0.43 






NGC 5471 


0.25 






NGC 6822 I 


0.02 


0.29 


0.68 


NGC 6822 III 


0.00 


0.32 


0.74 



Table 9. Differences in the determination of emission line fluxes 
between our photometry and values available in the literature. 



Aother = Fh^'c) ~ F H(p%) including available niters for 
each comparison. Hence, positive differences indicate our de- 
termined fluxes are stronger. 

6.2 Comparison with K79 

K79 obtained photographic plates (on Kitt Peak's 2.1 m and 
0.9 m telescopes) and photoelectric Ha photometry (Kitt 
Peak's 0.9 m and Manastash Ridge 0.8 m telescopes fitted 
with RCA C31034 GaAs phototubes) of a large number of 
Giant HII regions which include all of our sample. Plotted in 
Figure |l^ are our CCD photometry values on the x-axis vs. 
Kennicutt's on the y-axis. Some regions show good agree- 
ment while others are too far from the line that indicates 
identical values. The deviant regions are almost all from the 
same parent galaxy, M 101. Although these regions have an 
integrated flux that is within the range of the other regions, 
they have a much concentrated profile, with a high flux den- 
sity in the inner core as seen from the radial flux density 
profiles (figures |l| to MMTG argue that there could be 
some problem with saturation in Kennicutt's photometry, 
on the basis of a comparison with their H/3 photometry, and 
a low Ha to H/3 ratio. 

The differences in the Ha flux measured for regions 




a- 



H1.5 - A 




o 

-12.5 - 

i i i i l i i i i l i i i i l i 

-12 -11,5 -11 -10.5 

Emission line flux (This work) 

Figure 16. Comparison between our emission line photometry 
and previous ones available in the literature. Different symbols 
refer to different sources for comparison as follows: squares for 
K79, circles for M78, triangles for MMTG and the star symbol 
for Curchwell and Goss (1999). The dashed line shows the location 
of identical flux determinations in the plot. A shaded box at the 
lower right corner shows an approximate range of photometric 
errors. 

NGC 604 and NGC 595 are small, but are slightly larger 
for NGC 588 and NGC 592, all of them from the same host 
galaxy, M33. It is mentioned in K79 that photometry for 
lower brightness regions was obtained from photographic 
plates, which would enhance non linearity problems. 

Errors could also arise from any systematic but unde- 
tected difference between the sky conditions during our ob- 
serving run. This is a possibility, as most of the GEHRs in 
M101 were observed in the second half of the run. How- 
ever, the scatter observed for the standard stars in those 
nights is much smaller than the differences in flux, and al- 
though regions NGC 595 and NGC 5471 were observed on 
the same night, only NGC 5471 shows a discrepancy between 
the fluxes calculated in this work and those from Kennicutt. 



6.3 Comparison with M79 

There are five regions in common with the H/3 photometry 
performed by M79. In his work, the H/3 emission fluxes were 
determined with a two-channel photoelectric photometer at- 
tached to the 1.5m telescope at the Palomar Observatory. 
The flux calibration was performed through the observa- 
tion of three standard subdwarf F stars. Melnick estimates 
a conservative error of ±0.1 mag. although there could be 
somewhat larger systematic errors due to differences in an- 
gular sizes of the nebulae and its effect on the illumination 
of the photocatode. 

Figure |l(] shows a comparison between our H/3 fluxes 
and the ones found in his work. The five regions in common 
are plotted with open circles. As it can be seen from the 
figure, there is a systematic difference between both results, 
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being our values brighter by a factor of 2. A similar system- 
atic difference can be seen with the H/3 photometry which 
we will discuss below. 



6.4 Comparison with MMTG 

The comparison with MMTG photometry is more difficult, 
as there are five regions in common too, and there are no 
details about the photometric determinations in their pa- 
per. In figure jl6| we can see that the overall agreement is 
good, although the discrepancy for one region (NGC 588) is 
somewhat larger than expected. As it was already discussed 
in Section |^, our H/3 fluxes might be larger than expected 
when compared to Ha fluxes, as the ratio very rarely ex- 
ceeds the theoretical value 2.82, expected for no extinction. 
Furthermore NGC 588 shows the smallest ratio (2.34) which 
suggests an overestimation of our H/3 flux for that region. 

MMTG show in their work a very good agreement with 
M79 photometry, which might seem puzzling as our data 
agrees only with one of them. This arises from the fact that 
the agreement between the previous studies holds for the 
larger sample in common, but not for the subset of the re- 
gions included in our sample. 

Still, our H/3 photometry is not as reliable as expected. 
Keeping this in mind, we will base most of our further anal- 
ysis on the Ha photometry. 



6.5 The case of NGC 604 



Figurel7 



In a recent paper Churchwell & Goss (1999) study the ex- 
tinction within NGC 604 based on their own radio wave- 
length data from the Very Large Array (VLA) and CCD Ha 
photometry performed on an HST image from D. Hunter. 
We will further discuss their findings later in this paper. 
Their determined total Ha flux for the region, (3.9 ± 0.4) x 
10 -11 ergs cm~ 2 s" 1 is in agreement within the errors with 
our value of (3.1 ±0.4) x 10 -11 . 



7 EXTINCTION MAPS 

From the overall fluxes determined for each region we have 
calculated the ratio between the Ha and H/3 fluxes, in or- 
der to estimate the extinction affecting the measured Ha 
flux. However, it is quite well known for star forming regions 
that the nebular characteristics are far from homogeneous. 
This applies to the extinction properties of the interstellar 
medium, which combines neutral, ionised and molecular gas, 
plus dus t in d ifferent geometrical distributions. Churchwell 
& Goss (1999) have used the ratio between radio continuum 



and Ha to map the extinction properties in NGC 604 and 
analyse their correlation with the CO distribution from the 
work of Viallefond et al. (1992). A lo w spa tial resolution map 
was constructed by Mai'z Apellaniz (1999) using longslit flux 



calibrated spectra taken at different positions in NGC 604. 

In order to obtain a 'viewable' version of the flux ratio 
map, the calibrated nebular images had to be filtered to 
remove low signal pixels that introduce large scatter when 
performing the pixel to pixel division. For the sake of clarity, 
we have removed all pixels which had flux lower than seven 
times the background noise. After performing the division, 
the spatial variations were smoothed filtering the resultant 



Figure 17. Ratio of Ha to H/3 flux for NGC 604. The image has 
a larger scale than the nebular filter image, although it shares the 
coordinates origin 



Figurel8 



Figure 18. Same as figure [n] for NGC 595 



image. We chose a circular median filter with a four pixel 
radius, by means of the IRAF task rmedian. Although this 
introduces a degradation in the image spatial resolution, it 
helps to visualize the main trends in extinction changes. 

From our data, we can produce extinction maps 
for NGC 604, NGC 595, and both observed regions in 
NGC 6822. These are shown in figures |l| to 

The extinction varies strongly throughout most of the 
GEHRs (the Ha/H/3 flux ratio varies from 2.5 up to 7), 
which can be readily explained when considering the dra- 
matic changes in the distribution of obscuring material in 
these star forming regions. NGC 6822 I turns out to be 
a singular case, as the Ha/H/3 flux ratio remains constant 
and close to the theoretical expected value. This can be in- 
terpreted as a lack of obscuring matter (such as molecular 
clouds) which might suggest there is no material left over 
for further, or ongoing, star formation. We will discuss this 
issue in section B.4. 



8 THE LUMINOSITY - VELOCITY 

DISPERSION (L - cr) RELATION FOR 
GEHRS 

8.1 Luminosities 

With the observed fluxes obtained in the previous section 
it is possible to determine the luminosity by correcting for 
the distance d to the parent galaxy. The correction fac- 
tor is simply And 2 , with d expressed in cm. Distances are 
usually expressed in distance modulus (m — M) related by 
M — m + 5 — 5 log(D), where D is the distance in par- 
sees. The determination of accurate distance to the galaxies 
of the Local Group has been the central aim of several re- 
search projects. Distances to NGC 6822, M 33 and M 101 
have recently been determined with primary indicators, such 
as Ceph eids, by the HST Ext ra galactic Distance Scale Key 
Project (Freedman et al.1991), (Kelson et al.1996). We have 
used the cepheid distance scale for the galaxies hosting the 
GEHRS included in this study and listed in Table [t(]. This 
indicator also agrees well with other determinations such as 
the tip of the red giant branch and the planetary nebul a 
luminosity function, as discussed in Kennicutt et al. (1998). 
Note that the estimated distances to the host galaxies in- 
clude a considerable source of uncertainty which adds to the 
flux errors determined above. 

In order to compare with previous dete rmina tions, we 
have included results from M79, Hippelein (|l986|) (HH86), 



Figurel9 



Figure 19. Same as figure [| for NGC 6822 I 
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Figure20 



Figure 20. Same as figure [n] for NGC 6822 III 



Galaxy 


(m - M) PL 


log(47rd 2 ) 




&MMTG 


M 33 


24.63 ± 0.09 


49.93 


-0.23 


0.18 


M 101 


29.35 ± 0.17 


51.82 


+0.26 


0.06 


NGC 6822 


23.62 ± 0.20 


49.53 


+0.09 


0.13 



Table 10. Distance to the galaxies hosting the GEHRs included 
in this study. 



and MMTG. However, it must be kept in mind that each 
study has made use of the available distances from different 
sources and discrepancies in these distances introduce differ- 
ent corrections. We have included in Table [n] the correction 
to the logarithmic luminosity calculated with the distances 
adopted in this work, and the differences with the same fac- 
tors derived for the distances used in previous works. It can 
be readily seen that the distance to the galaxy is a consid- 
erable source for disagreement in luminosities. The magni- 
tude of these disagreements on the distance scales used and 
the problems that may have affected previous photometric 
studies, as described in the previous section make any com- 
parison of the calculated luminosities almost impossible. 



8.2 Velocity dispersions 

Velocity dispersions have been independently obtained by 
several groups analysing the kinematics of GEHRs. These 
efforts can be split according to the spatial resolution of the 
instruments used for the estimations into large and small 
scale observations. I will discuss in this section the large 
scale observations, as these are directly linked to the over- 
all luminosity measured for the HII regions. Melnick (1977), 



measured line profiles and determined turbulent velocities 
for 



a sej l of GEHRij 

interfer DTTTu (XT 



rtr 



in nearby galaxies 
Paluniai 



using a Fabry-Perol 
bsei vaLuiy — when investigating 

the relation of the velocity dispersion with the linear diam- 
eter of giant HII regions. The turbulent velocity a is deter- 
mined from the observed line profile a Q after correcting from 





M79 


HH86 


MMTG 


This work 


GEHR id. 


Ih/3 


I Ha 


Ih/3 


I Ha 


NGC 588 


38.08 


38.61 


38.25 


38.47 


NGC 592 




38.42 


38.44 


38.33 


NGC 595 


38.64 


39.18 


38.55 


38.95 


NGC 604 


38.83 


39.64 


38.80 


39.42 


NGC 5447 




39.82 


39.59 


40.45 


NGC 5461 


39.93 


40.12 


39.77 


40.62 


NGC 5462 


39.72 


39.80 


39.65 


40.43 


NGC 5471 


39.73 


40.17 


39.85 


40.62 


NGC 68221 


38.46 


38.36 


37.59 


38.21 


NGC 6822III 


38.30 


38.44 


37.59 


38.30 



Table 11. Luminosities available in the literature, plus those 
calculated in this paper. 



GEHR id. 


HH86 


RAJ86 


MMTG 


NGC 588 


12.9 ± 0.4 


14.9 (21.1 ± 1.0) 


12.8 


NGC 592 


13.6 ± 1.0 


12.1 (17.1 ± 1.0) 


12.6 


NGC 595 


17.8 ± 0.8 


19.2 (27.1 ± 1.6) 


17.2 


NGC 604 


17.8 ± 0.3 


16.3 (23.1 ± 0.8) 


16.8 


NGC 5447 


20.3 ± 0.6 


21.2 (30.0 ± 1.2) 


19.3 


NGC 5461 


21.5 ± 0.3 


24.7 (35.0 ± 1.3) 


20.3 


NGC 5462 


24.2 ± 0.4 




21.8 


NGC 5471 


22.4 ± 0.4 


21.3 (30.2 ± 1.0) 


21.0 


NGC 68221 


10.5 ± 0.2 


9.5 (13.4 ± 0.8) 


9.8 


NGC 6822III 


10.3 ± 0.3 


10.4 (14.7 ± 0.9) 


9.5 



Table 12. Velocity dispersions from the literature, of the ionised 
gas of the GEHRs included in this study. We have listed the values 
of the a dispersions published, with the exception of the values 
for RAJ86, which we scaled from the beta width listed in their 
paper, included between parentheses in column 3. 



the instrumental profile Oi and the thermal contribution to 
the observed width of the emission line <jt- Hence: 



2 2 2 2 

a = a a - <Tj - a T 



(16) 



HH86 obtained newer F-P interferometric data, including 
a larger sample of region s with dissimilar results. Indepen- 
dently Roy et al. ( 198£ ) (RAJ86) obtained another set of 
velocity dispersions for extragalactic HII regions. The agree- 
ment of both sets of measurements for the regions included 
in our study can be easily seen in Table ha. In a later 
analysis MMTG combined the data available for Ha and 
[OIII] from HH86 and characterised the gas kinematics with 

CHa + C[OIII] 

a= 2 -■ 

There are several other studies of spatial variations 
in the kinematics of these large regions. These were done 



mostly with echelle spectr oscopy (Skillman & Balick (1984), 
Chu & Kennicutt Jl994| )) and with the TAURUS Fabry- 
Perot Imaging Spectrograph on the 4 m. William Herschel 



Telescope in L a Palma (Muno z-Tun on et al. ( |1995; MGC95|), 
(h.996; MTCTl ), Yang et al. ( |199^ )). TAURUS~ paylor and 
Athertonl980D allows to obtain simultaneously high resolu- 



tion spectroscopic information with its etalon grating and 
seeing limited imaging. These studies added new informa- 
tion regarding the smaller scale kinematics of these regions, 
showing there is a combination of shell structures, loops and 
supersonic motions down to the smallest observa ble d etail. 
We will discuss these important results in section S.4. 



8.3 The \og(L(Ha)) vs log(cr) regression 

As already mentioned, several studies analysed the existence 
of a correlation between the total luminosity emitted by a 
giant HII region in one of its Balmer lines and the velocity 
dispersion of the gas measured from the width of such line. 
The slope derived for this regression is a crucial parameter 
to find an explanation for the source of such a relation. In 
this section we will analyse and discuss the existence of this 
correlation on the basis of our photometry. 

We have combined the luminosities obtained in the pre- 
vious section from our CCD photometry and the kinematical 
data used in MMTG in order to obtain an improved set of 
data for analysing their relation. The existence of a corre- 
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1.1 1.2 
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Figure 21. Plot of the logarithm of the Ho luminosity against 
the logarithmic velocity dispersion. A linear fit to the data is also 
drawn. 



lation was analysed by performing a least squares fit to the 
available data, accounting for errors in both coordinates. Er- 
rors were estimated to be about 10% in the kinematics and 
20-25% in the luminosities. In the logarithmic plane, a linear 
regression is found, where 

\og(F(Ha)) = 29.6(±0.7) + 8.2(±0.6) x log(a). (17) 

The data is plotted in Figure |2l] together with the line 
showing the fit to the points. From the plot it can be readily 
seen that although there is a correlation, there is also a large 
scatter in the data. Deviations from the fit are too large for 
the expected uncertainties (a \ 2 test on the fit yields a prob- 
ability almost zero) and large enough to seriously limit its 
use to derive accurate luminosities from measured emission 
line widths. 

It seems that, under the evidence of the new photometry 
here presented, the uncertainties in previous determinations 
of giant HII region luminosities were larger than estimated. 
Also, there is increasing evidence that there is an age spread 
between GEHRs and that this could modify noticeably the 
observed global parameters. However, our knowledge of the 
finer details of some of these regions can help us select a 
more homogeneous sample. 



8.4 The young GEHRs 

The small scale kinematic studies brought new information 
to the analysis of GEHRs. Most of them showed that the 
supersonic profile obtained with global spectroscopic data is 
formed by the contribution of different regions within each 
GEHR with different kinematic behaviour. It was also found 
that this mixture of different contributions was not uniform 
for all HII regions; some regions have most of their flux com- 
ing from supersonic cores, others have a larger contribution 
from bubbles and shells. We have used this information to 
select a more homogeneous set of GEHRs out of our sample. 



Particularly for the farthest regions, for which the angular 
resolution makes any guessing of its structure impossible. 
Based on the assumption that the large shells and bubbles 
are formed at later stages of the evolution of these star form- 
ing regions is that we will refer to those regions that do not 
show evidence of loops and filaments as the young GEHRs. 
We will use spatially resolved kinematical data to disentan- 
gle the multiple blobs that form the M 101 GEHRs and 
identify those who show supersonic velocity widths. Includ- 
ing the available information, and making use of the ex- 
cellent seeing conditions for our photometric data we have 
re-analysed all the regions in our JKT photometry sample. 

8.4.1 M 33 GEHRs 



The galax y M 33 is situated at about 840 kpc (Freedman 
et aJ.1991 ). The whole emission line region of the 30 Do- 
radus nebula in the Large Magellanic Cloud, which is the 
one for which we have the best detailed information, would 
occupy only 1 arcmin, instead of the 15 arcmin it fills at 
its actual distance. However, at this distance it is still pos- 
sible to resolve the presence of large bubbles, shells or any 
other evidence of nebular disruption and combine it with 
other information regarding its age or evolutionary stage of 
its components. 

NGC 588 

This region shows clear evidence, as seen in figure ^ that 
there is a compact core, but it also has at least two strong 
ring-like structures that contribute to the emitted flux. Kine- 
matic data available from MTCT indicate the presence of 
shells with intensity comparable to the rest of the kinematic 
core. From the expected traces of evolution of GEHRs we 
consider this region as an older region which has been swept 
by strong winds and supernova explosions. 

NGC 592 

The Ha image shows two nuclei dominating the emis- 
sion of the core of the region. A filamentary nebula surround- 
ing the core and extending out to large distances, makes it 
difficult to find an asymptotic value to determine the total 
flux (see Figure ^|). There is no data about the small scale 
kinematics of this region, although its structure reveals a 
strong interaction with the external environment, suggest- 
ing an older age too. The presence of two condensed knots 
makes it an interesting target for analysing their kinematical 
properties with TAURUS-like Fabry-Perot interferometers. 

NGC 595 

This region has a structure quite similar to that of 
NGC 592. There are a number of dispersed strong knots, 
visible in figure ^| There is no data available from small 
s cale kinematics for this region. However, Malumuth et al. 
have obtained multiband HST/WFPC2 images of the 



ionizing stellar cluster. They derive an age of 4.5 Myr which 
indicates it is a somewhat evolved cluster consistent with 



the presence of Wolf-Rayet stars. (Drissen et al.1993) 



NGC 604 

The region NGC 604 shows an extended structure, and 
some small loops can be identified within the nebula. There 
is, however no strong evidence of disruption of the nebula, as 
the stellar cluster remains enshrouded by the nebula in the 
Ha image shown in figure ^. A comprehensive investigation 
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Figure22 



Figure23 



Figure 22. Expanded view of the components of NGC 5447. The 
most important knots are marked and labelled. The circle drawn 
illustrates the apparent size that the 30 Dor nebula would have 
at the distance of M 101. 



of t he kin ematics of NGC 604 was carried out by Yang et 
al. (1996). In their work they gathered longslit echelle spec- 
troscopy and TAURUS data and detected expanding shells 
together with an overall presence of broad profiles in ev- 
ery position within the region. They plotted the measured 
FWHM of the profiles versus their intensity and found that 
the brightest profiles define a narrow band in such diagram. 



This band, with a median value of ~ 36 km s include 
most profiles of the distribution, and dominates the total 
observed flux. TAURUS data in MTCT shows a similar dis- 
tribution in the / vs. a diagram. It was proposed in both 
papers that gravitation would be responsible for the pro- 
files in the narrow band, and hence dominat e the kinematics 
of the region. Gonzalez Delgado & Perez ( 2000 ) combined 



archival HST data, together with spectra from the Interna- 
tional Ultraviolet Explorer (IUE) and optical spectra from 
the William Herschel Telescope (WHT) to analyse the evo- 
lutionary state of the region. They combined the available 
information from stellar wind resonance lines, nebular op- 
tical emission lines and Hel absorption lines present in the 
spectra of NGC 604 to derive a 3 Myr age for a single burst 
that originated the central ionising cluster. All these facts 
indicate that NGC 604 can be considered as a young GEHR. 



8.4.2 M 101 GEHRs 



The spiral gala xy M 101 lies at about 7.4 Mpc from us (Kel- 
son et al.1996). This is ten times more distant than M 33 



and consequently with a linear resolution one order of mag- 
nitude smaller. The apparent size of the 30 Dor nebula at 
this distance would be of only 6 arcsecs. It is very difficult to 
base an analysis simply on the distribution of the emission 
flux, but we will rely on the information available from other 
kinematic studies, whenever available. 

NGC 5447 

This region, as almost all regions in M 101 has already 
been classified as ' multiple' due to the prese nce of several 
bright components ( Lawrie and Kwitterl982 ). We have in- 
cluded an expanded view of this region in figure [22] and la- 
belled the most conspicuous of these components or 'knots'. 
Of particular interest is knot A, which is the strongest of 
all. There is no detailed kinematics for this region, although 
HH86 took spectra of the northern and southern part and 
found them to be almost identical. Considering this informa- 
tion we performed aperture photometry, taking increasing 
apertures centred in NGC 5447 A, using the local diffuse 
nebula to estimate the background level. In this way the 
net Ha luminosity of this particular knot is calculated, and 
found to be ~ 7.5 x 10 -13 ergs s _1 . 

NGC 5461 

The multiplicity of NG C 5461 had already been sug- 
gested by Kennicutt (1984) and analysed by MGC95 with 
TAURUS. The knots, as identified by MGC95 are marked on 
our Ha image in figure E3l It can be predicted from the im- 



Figure 23. Same as fig^ for NGC 5461. T he nom enclature used 
for the knots follows that of Munoz-Tunon (1994). 



Figure24 



Figure 24. Same as fig|22j for NGC 5462. 

age that knot A dominates the overall flux. TAURUS data 
revealed that 'well-behaved' Gaussian emission lines show 
a velocity dispersion of 25 km s _1 . Furthermore, these su- 
personic sigma values are restricted to knots A and B only. 
Knots C and E are kinematically detached from the region. 
This can be guessed from the figure, although the definitive 
confirmation is given by the resolution available with 2D 
spectroscopy. 

In a similar way as described for NGC 5447 A we have 
calculated the Ha luminosity emitted by NGC 5461 A and 
found it to be ~ 2.6 x 10~ 12 ergs s" 1 . 

NGC 5462 

The complex structure of this giant HII region is shown 
in figure |24|. Two major knots, labelled A and B are much 
brighter than the remaining knots, but this particular region 
does not show a single dominant feature as is the case of the 
previous ones. There is no information available regarding 
the kinematical behaviour of these knots as to assign the 
characteristics of the global profile to any of these knots in 
particular. The approximate luminosities of knots A and B 
emitted in Ha are 3.5 and 4.5 xlO -13 ergs s _1 respectively. 

NGC 5471 

NGC 5471 also shows multiple components. Skillman 
analysed spatial variations for this region and iden- 

In that figure 
detected 



( |L985| ) 



tified 5 knots which are labelled in figure 
it can also be seen another weak component, A' 
under better seeing conditions by MGC95. This latter study 
suggests that a large fraction (80%) of the region could be 
characterized by a constant emission linewidth of a ~ 20 km 
s" 1 . We have estimated the Ha luminosity emitted by knot 



A to be ~ 10" 



ergs s 



8.4.3 NGC 6822 GEHRs 

NGC 68 22 is an irregular gala xy, situated at approximately 
530kpc ( Kennicutt et a!.199£ ). It is the closest one from the 
sample studied here, and therefore the one that allows the 
best linear resolu tion, at 2.5 pc arcsec - . O'Dell, Hodge and 
Kennicutt (1999) obtained HST photometry (using mainly 
emission line filters) and determined for both stellar clusters 
associated with the two regions analysed in this paper, an 
approximate age of 4 Myr. 

NGC 6822 I 



according to Skillman 





Figure25 
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GEHR name 


Age 


Kincmatical core 


NGC 588 


> 4.5 Myr 




NGC 592 


> 4.5 Myr 




NGC 595 


4.5 Myr 




NGC 604 


3 Myr 




NGC 6822 I 


4 Myr 




NGC 6822 III 


4 Myr 




NGC 5447 




NGC 5447A 


NGC 5461 




NGC 5461A 


NGC 5462 




nil 


NGC 5471 




NGC 5471A 



Table 13. Summary of analysis performed over the GEHRs sam- 
ple to detect regions with signs of evolution or multiplicity. Ages 
are as derived from studies of the stellar clusters associated with 
the regions, or from an assumed evolution of the nebular fea- 
tures (see text). The identification of kinematical cores was made 
from spatially resolved kinematical data of the multiple regions 
in M 101. 



This region (figure shows a spherical concentrated 
bulk, and an extended loop (or pair of loops). There is no 
detailed kinematical data, but one can expect that the global 
emission profile is dominated by the central core. However, 
the extinction map reveals that there is almost no dust left. 
This behaviour is not observed in other regions that show 
recent star formation activity and might suggest some evo- 
lution effects, though not as strong as in other GEHRs. 

NGC 6822 III 

This region has the emission highly concentrated too, 
with a bright nucleus and a nebular envelope, with no evi- 
dence of shell-like features (see figure ^) . 



8.4-4 The \og(L(Ha)) vs. log(cr) for the young and single 
GEHRs 

The analysis indicates that the GEHRs included in this 
study do not necessarily share the same global properties. 
This is evident in the case of total luminosity and global ve- 
locity dispersion. From the expected evolution of the kine- 
matics of gi ant HI I regions in the frame presented by Tenorio 
Tagle et al. ( |l993| ) and MTCT young GEHRs (those that do 
not show signs of evolution) should be governed by the grav- 
itational supersonic motion induced by the stirring stars. To 
test this hypothesis, we have selected a sub sample that in- 
cludes only those regions which do not show evident signs of 
evolution, such as loops or shells. We have also included the 
dominant cores of the multiple regions in M 101, as their 
sizes are comparable to the 30 Dor Nebula, suggesting by 
analogy, that these cores are the most prominent members 
of an agglomeration of HII regions. 

The selection process is summarised in Table jl3| Al- 
though the age determinations are very dependent on the 
used evolutionary tracks for each case, we are discarding 
regions that are 4.5 Myr old, or seem to be even older. Re- 
garding the multiplicity, we are only including knots that 
have properties that suggest the existence of a kinematical 
core. 

No matter what quantitative criteria is established, the 
final decision turns out to be subjective at this stage. There 




1.1 1.2 

log(velocity dispersion} 

Figure 26. The logarithm of the Ho luminosity against the log- 
arithmic velocity dispersion for the subsample described in this 
section. The regression line with slope 4.73 is also shown. 



is a compromise between discarding any region that shows 
signs of evolution and keeping a fair number of GEHRs as 
to attempt a regression analysis. Therefore, the short list 
includes: NGC 6822 I, NGC 6822 III, NGC 604, NGC 5471 
A, NGC 5461 A, and NGC 5447 A. A similar plot as the 
one presented in figure |2lj but for the mentioned subsample 
is shown in figure Eq. In it we have included the linear fit to 
the points, which was found to be 



log{F(Ha)) = 33.6(±0.6) + 4.7(±0.9) x log(a) 



(18) 



The slope and zeropoint are - within the errors - in agree- 
ment with the values found by Fuentes-Masip et al. (200C) 
for their sample of high surface brightness giant HII regions 
in NGC 4449, supporting a virialised system as the origin of 
the supersonic motions. 



9 DISCUSSION 

This paper presents the results of new Ha and Hj3 CCD 
photometry of ten Giant HII Regions lying in Local Group 
galaxies. The comparison of these fluxes with previous work 
suggests some systematic errors affecting the flux determi- 
nations made with photoelectric detectors. 

A serious problem is the lack of a uniform method for 
flux calibration and for that we have introduced the use of 
spectrophotometric standards in order to rigorously calcu- 
late the continuum flux correction. Accounting for the fact 
that the filter used is broader than the emission feature is 
very important, and it is generally not discussed in papers 
dealing with this issue. Accurate knowledge of the emission 
line characteristics are crucial in order to calculate its po- 
sition within the filter transmission curve and its relative 
width. All these make the narrowband photometry of emis- 
sion line objects a very difficult one, specially when they 
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are extended objects and aperture spectrophotometry is not 
possible either. 

The analysis presented here has two important although 
apparently dissimilar results: 

• Giant Extragalactic HII Regions are far too heteroge- 
neous objects as to blindly include them in a sample for use 
as an accurate (rms < 0.2 \og(L(Ha)) distance indicator. 

• Young GEHRs and the kinematical cores of GEHRs 
included in this study do follow a tight relation in the L(Ha) 
vs. a plane. 

The first item does not preclude the use of GEHRs as 
a distance indicator, but establishes that a careful selection 
based on the kinematical properties - or ages - of these is 
a must, as the luminosities of the young GEHRs and kine- 
matical cores seem to correlate very well with the veloc- 
ity dispersion of the region. The use of the total Ha (or 
H/3) luminosity of the region has not eliminated the prob- 
lem of dealing with other parameters of the nebula, as it 
was hoped it would. What were thought to be 'supergiant' 
GEHRs (such as NGC 5461 and NGC 5471) should be bet- 
ter considered as a group of HII regions, out of which one 
(or several) of them might be a true GEHR. The emission 
line flux originated in non-giant HII regions can be partially 
responsible for the observed scatter in the regression. 

The second item mentioned above indicates that an- 
other factor that turns out to be very important is the age, or 
evolutionary stage, of the regions. It has already been show n 
for the particular case of 30 Doradus (Selman et al.1999), 
that older bursts of star formation can coexist in an active 
star forming region. Hence, the presence of clear signatures 
of recent massive star formation do not necessarily indicate 
that we are looking at an extremely young HII region. Al- 
though the most recent burst of stars will be responsible 
for the bulk of the ionisation of the gas - and hence of the 
emitted flux - the effect of previous stellar generations could 
have modified substantially the interstellar gas kinematics. 

New observational efforts should be directed to the use 
of imaging spectrographs to analyse Giant HII Regions in 
the Local Group to learn more about the kinematical char- 
acteristics of the emitted flux. Other global properties from 
HII regions, such as the equivalent width Wh/3 of the H/3 lin e 
known to vary with the evolutionary stage (Dottori (1981), 
Copetti et al. (1986)) can be used to estimate the age of the 
GEHRs. A calibration of a larger sample will be decisive to 
prove what is suggested from the data now available and to 
provide a link with what is observed for starburst galaxies. 
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